A06 35451 


Attribute  Sampling  Plans  Based  on  Prior  Distributions  and  Costs. 


e 

By 

A.  Ha Id 

University  of  Copenhagen. 

^1.  A  model  with  linear  costs. 

Let  N  and  n  denote  lot  size  and  sample  size,  and  let  X  and  x  denote  number  of 
.defectives  in  the  lot  and  the  sample,  respectively.  The  acceptance  number  is  denoted 

1 

by  c. 

Let  the  costs  be 


-  D 

'Hi 
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and 


nSj+  xS^l-  (N-njA^  +  (X-x)A2 

for  x  £  c 

(1) 

nSL+  xS2+  (N-n)R1  +  (X-x)R2 

for  x  >  c. 

(2) 

1 


The  interpretation  of  the  six  cost  parameters  depends  on  the  kind  of  inspection 
envisaged,  i.e.  whether  inspection  is  a  consumer's  receiving  inspection,  a 
producer's  inspection  of  finished  goods,  or  "internal  inspection"  by  delivery 
of  goods  from  one  department  to  another  within  the  same  firm.  The  cost  parameters 
may  have  quite  different  values  when  considered  exclusively  from  a  producer's  or 
a  consumer's  point  of  view  because  certain  costs  are  borne  primarily  by  one  of  the 
parties  involved.  The  values  of  the  cost  parameters  also  depend  on  whether  the 
inspection  is  rectifying  or  non-rectifying,  destructive  or  non -destructive.  In  the 
following  the  two  cost  expressions  are  discussed  and  a  few  examples  of  interpretation 
are  given. 

Costs  associated  with  the  sample,  nS^+  xS2,  for  brevity  called  "costs  of  sampling 
inspection",  consist  of  two  parts:  one  part,  nS^,  proportional  to  the  number  of 
items  in  the  sample  so  that  includes  sampling  and  testing  costs  per  item,  and 
another  part,  xS2,  proportional  to  the  number  of  defectives  in  the  sample,  i.e. 

S2  denotes  additional  costs  for  an  inspected  defective  item.  If  defective  items 
found  in  the  sample  are  repaired,  say,  then  S2  includes  repair  costs  per  item. 

"Costs  of  acceptance"  are  similarly  composed  of  a  part,  (N-n)A.,,  proportional  to 
the  number  of  items  in  the  remainder  of  the  lot,  and  another  part,  (X-x)A2, 
proportional  to  the  number  of  defective  items  accepted.  Whereas  A^  usually  will 
be  zero  or  negligible,  A9  will  often  be  considerable.  If  accepted  items  are  used 
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parts  in  an  assemply  operation,  say,  may  include  the  manufacturing  costs 
(or  the  price)  of  an  item,  the  costs  of  handling  the  defective  item  in  assembling 
and  disassembling,  and  the  damage  done  to  other  parts  used  in  the  assembly.  In 
case  of  inspection  of  finished  goods  may  include  costs  of  repair,  service  and 
guarantees  plus  loss  of  good-will, 

"Costs  of  re  .lection"  consist  of  a  part,  (N-n)R^.  proportional  to  the  number  of 
items  in  the  remainder  of  uhe  lot,  and  another  part,  (X-x)R2,  proportional  to 
the  number  of  defective  items  rejected.  Re  lection  is  here  taken  in  a  broad  sense 
meaning  only  that  the  lot  cannot  be  accepted  according  to  the  sampling  plan  used. 
Rejection  may  therefore  lead  to  sorting,  price  reduction,  scrapping,  or  salvaging. 

If  rejection  means  sorting,  say,  then  R^  includes  sorting  costs  per  item  and  R^ 
denotes  additional  costs  for  defective  items  found,  for  example  costs  of  repair 
or  replacement. 

It  is  obvious  that  from  a  practical  point  of  view  it  will  in  general  be  easiest 
to  obtain  Information  on  the  values  of  the  cost  parameters  in  the  case  of  "internal 
inspection". 

Let  £^00  denote  the  (prior)  distribution  of  X,  i.e.  the  distribution  of  lot  quality. 
It  is  assumed  that  this  distribution  is  a  mixed  binomial  distribution,  i.e. 

£„<X)  -  }$pVXdW<p)  (3) 

where  W(p)  denotes  a  cumulative  distribution  function  (independent  of  N) . 

Drawing  a  sample  without  replacement  from  each  lot  (hyper geometric  sampling)  and 
computing  the  average  costs  we  find 

1 

K(N,n,c)  -  /  K(N,n,c,p)dW(p)  (4) 

0 

where 

K(N,n,c,p)  =  n(Sj+  S2p)  +  (N~n) ((Aj+  A2p)P(p)  +  (Rj+  R2p)Q(p)),  (5) 

and  c 

P(p)  =  B(c,n,p)  =  E  (£)pXqn‘X,  Q(p)  -  l-P(p).  (6) 

x=0 

For  convenience  the  frequency  function  corresponding  to  W(p)  will  be  called  the 
distribution  of  the  process  average  or  the  distribution  of  p  as  distinct  from 
f^(X)  which  gives  the  distribution  of  X/N,  i.e.  the  distribution  of  lot  quality. 

Starting  from  (5)  we  introduce  the  three  cost  functions 

k8<p>  "  sx+  S2P>  kfl(P)  *  A!+  a2P>  kr<P>  *  Ri+  V' 
the  (economic)  break-even  quality. 

Pr=  (Rx-  A1)/(A2-R2),  0  <  pr  <  1,  and 

k  (p)  ■  k  (p)  for  p  £  p^  and  k  (p)  for  p>  p  . 
m  a  r  r  r 


(8) 


The  function  k  (p)  gives  the  unavoidable  (minimum)  costs,  i.e.  the  coots  correspond* 

m  of  informatioi 

ing  to  the  situation  where  perfect  knowledge  of  quality  exists  without  costs  and 

all  lots  are  classified  correctly  on  the  basis  of  the  corresponding  process 

average,  viz.  accepted  for  pS  pf  and  rejected  for  p>  p^. 

Averages  over  the  prior  distribution  are  denoted  by  k  ,k  ,  etc.,  and  costs  per 

^  8  Si 

item  are  denoted  by  k,  and  costs  per  lot  by  the  corresponding  K,  i.e.  K  *  Nk. 

The  average  costs  for  the  three  cases  without  sampling  inspection,  i.e.  the  cases 
where  (a)  all  lots  are  classified  correctly,  (b)  all  lots  are  accepted,  and  (c) 
all  lots  are  rejected,  then  become  k^k  ,  an<*  respectively.  These  cases  are 
useful  "reference  cases"  since  sampling  inspection  is  justified  only  if 
k-k  <  min{k  -  k  ,k  -  k  },  where  k  *  K(N,n,c)/N. 

Case  (a)  will  usually  be  considered  as  the  basic  reference  case  and  average  costs 
for  other  cases  will  therefore  be  reduced  by  k^,  since  k^  represents  the  average 
fixed  costs  per  item  which  will  be  incurred  irrespective  of  the  decision  made. 

The  cost  differences 

1  pr 

k  -  k  -  /  (k  (p)-k  (p»dW(p)  and  k-k  «  /  (k  (p)-k  (p))dW(p) 
am^a  r  r  m  q  r  a 

represent  average  decision  losses  in  case  (b)  and  (c)  respectively,  and  k  -  k 

8  IQ 

represents  the  average  "loss"  by  inspection. 

From  (4)  and  (5)  we  find 

Pr  1 

K-K  «  n(k  -k  )+(N-n)(A.-R  )(/(p  -p)Q(p)dW(p)  +  f(p-p  )P(p)dW(p) ),  (9) 

m  sm  <  <  .  c  r 

O  p 

*r 

the  two  terms  giving  the  average  costs  of  sampling  inspection  and  the  average 
decision  losses,  respectively. 


Putting  R(N,n,c)  *»  (K(N,n,c)-K  }/(k  -  k  ),  we  find 

Pr  1 

R  -  n  +  /( p  -p)Q(p)dW(p)  +  /( p-p  )P(p)dW(p)),  (10) 

*s  *m  0  p 

r 

the  two  terms  again  giving  the  costs  of  sampling  inspection  and  the  average 
decision  losses,  respectively,  but  here  using  the  average  costs  of  sampling 
inspection  (minus  k^)  per  item  in  the  sample  as  economic  unit. 

2.  Results  for  double  binomial  prior  distributions. 

The  simplest  theory  of  sampling  inspection  is  based  on  the  assumption  that 
lots  submitted  for  inspection  originate  from  one  of  two  possible  quality  levels, 
p^  and  p p^  <  p^>  and  that  the  process  average  p^^  occurs  with  probability  w^, 
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consequently  with  probability  w2«  1-w^.  The  standardized  cost  function  (10) 

1  then  takes  the  form 

R  =  n  +  (N-nXy^Cp^  +  72p(P2))  C11) 

^hcre  7j=  |Pi“PrlWi/(ps'Pra)  =  IW “MV  K/(W *  1  *  l’2’  and  Pm^lV’VV 
i.e.  R  depends  on  four  parameters  only,  viz. 

Since  P(p)  »  B(c,n,p)  denotes  the  operating  characteristic,  Q(p^)  and  PCp^)  equal 
the  producer's  and  the  consumer's  risks,  respectively. 

We  shall  discuss  ten  systems  of  sampling  plans  defined  as  follows: 

(1).  Bayesian  plans,  i.e.  plans  minimizing  R. 

Restricted  Bayesian  plans, i.e.  plana  minimizing  R  under  some  suitably  chosen 
restriction  on  the  operating  characteristic,  viz. 

)  (2).  Min  R  for  Q^)  -  a  or  P(p2)  =  (3.  p  q 

(3) .  Min  R  for  P(p_)  ■  1/2  where  p  =  (log  —)/(log  --2— -) . 

0  0  <12 

(4) .  Min  R  for  Q(p1>  *  a/N  or  P(P2)  =  P/N. 

(5) .  Min  R  for  P(p2>/Q(p1>  -  p. 

Plans  defined  by  two  risks,  vis* 

(6) .  Q(P1)  =  a/N  and  P(p2)  *  p/N. 

(7) .  Q(px)  “  «  and  P(p2)  -  0/N  (or  P(p2>  *  p  and  Q(px)  =  a/Nj. 

(8)  .  P(pQ)  «  1/2  and  Q(p1)  »  a/N  (or  P(p2>  -  P/N) . 

(9) .  Q(PX)  c  a  and  P(p2>  «  p. 

Finally  we  consider  percentage  inspection  defined  as 

(10) .  n  =  (iN  and  c  =  pQn. 

*  In  all  these  definitions  a, |3,p,  and  n  represent  suitably  chosen  positive  constants 
which  may  be  different  from  case  to  case. 

For  each  system  of  sampling  plans  it  has  been  shown  in  [lj  how  the  exact  solution 
may  be  obtained  and,  3ince  this  solution  is  an  implicit  one,  an  explicit  solution 
is  given  as  an  asymptotic  expansion  for  N  —>  oo  .  Tables  have  been  provided  in  [?.] 
and  [3]  . 

In  this  section  we  shall  assume  that  the  quality  distribution  of  submitted  lots  is 
a  double  binomial  and  compare  the  costs  of  the  various  systems  under  this  assumption. 
As  a  measure  of  efficiency  we  use  the  ratio  e(N,n,c)  =  R^(N)/R(N,n,c),  where  Rg(N) 
denotes  the  costs  of  the  optimum  (Bayesian)  plan  and  R(N,n, c)  denotes  the  costs 
of  the  plan  in  question. 

i 

(  The  systems  defined  by  (1)  -  (8)  fall  into  two  classes  depending  on  whether  both 

-1  ' 

risks  are  0(N  )  or  one  of  the  risks  is  constant  and  the  other  is  0(N  )  . 

I 
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The  first  class  contains  systems  (1),  (3),  (4),  (5),  (6),  and  (8).  Asymptotically 
the  relation  between  acceptance  number  and  sample  size  has  the  form 
c  ■  p^n  +  aj+  a^n  0(n  ^),  and  In  N  ■  cp^n  +  jin  n  +  *2n  0(n  which 

by  inversion  determines  n  as  a  function  of  N,  and  R  *  n  +  Bj+Bjn  *+  0(n  ),  the 
constants  pQ  and  <p^  being  the  same  in  all  cases  (depending  on  (p^pg)  only), 
whereas  the  remaining  constants  are  found  as  functions  of  the  parameters  in  the 
model  and  the  restriction.  This  means  that  n  =  0(ln  N)  and  that  the  average  decision 
loss,  (N-n)0'jQ(p1)+72P(p2))>  teRds  to  a  constant  (because  Q(pj)  and  P^)  are  0(N  1)). 


The  second  class  consists  of  systems  (2)  and  (7).  Asymptotically  the  relation 

4  l-i/2  -3/2 

between  acceptance  number  and  sample  size  has  the  form  c  ■  p.n  +  £  a,n  +  0(n  ), 

J  i»l  1 

Pj  representing  the  quality  level  having  a  constant  risk, 

1  •*  _i#p 

In  N  *  ip  n  +  r-ln  n  +  2  K.n1  '  ‘  +  0(n  ) ,  which  determines  n  as  a  function  of  N, 

■J  *•  4-1  1 

*•  ^  •1/2 

and  R  -  BN  +  (l-8)n  -{-  5^+  0(n"  7  ),  the  constants  cp^  and  6  being  the  same  in  the 
two  cases.  Because  of  the  constant  risk  all  relations  are  considerably  more  complicated 
than  for  the  first  class  and  R  becomes  0(N)  instead  of  0(ln  N) .  For  large  lots  it  is 
therefore  not  advisable  to  use  a  system  with  a  fixed  comsumer's  or  producer's  risk 
and  correspondingly  high  costs  as  compared  with  a  system  having  decreasing  risk'*. 


The  system  with  both  risks  fixed  and  the  system  with  percentage  inspection  both 
lead  to  R  ■  O(N)  and  asymptotically  they  have  the  same  costs  for  [a  ■  y^a  -|-  y^P. 

The  system  with  fixed  rioks  uses  a  fixed  sample  size  so  that  the  decision  los3  becomes 
of  order  N,  whereas  percentage  inspection  has  n  =  O(N)  and  a  decision  loss  of  order 


The  restricted  Bayesian  plans  with  both  risks  decreasing  and  the  corresponding  plans 
based  on  two  risks  have  an  economic  efficiency  tending  to  1  for  N  ->  co  as  compared 
to  the  Bayesian  plans.  (The  efficiency  pf  plans  having  at  least  one  risk  fixed  tends 
to  zero).  Tuis  result  means  that  wrong  values  of  the  weights  of  the  prior  distribution 
and  wrong  values  of  the  cost  parameters  have  a  secondary  influence  on  the  efficiency 

‘X*  * 

which  tends  to  1  if  only  (p^^)  are  correct.  If  wrong  values  of  (p^^)*  (PpP2)  say, 
are  used  for  finding  the  plans  then  the  efficiency  tends  to  e,  0  <  e  <  1,  if  and  only 

¥r 

if  Pj<  p^  <  p£<  p2>  otherwise  the  efficiency  tends  to  0. 

As  an  example  consider  a  case  with  p^=  0.01,  w^  0.85,  and  p2=  0.05,  W2S  0.15.  Let 
the  costs  of  sampling  inspection  be  0.40  (economic  units)  per  item  in  the  sample,  i.e. 
Sj«  0.40  and  82“  0,  the  costs  of  rejection  0.30  per  item  in  the  remainder,  i.e. 

Rj=  0.30  and  R^®  0,  and  the  costs  of  acceptance  per  defective  item  10.00,  i.e. 

Aj=  0  and  10.00.  It  follows  that  pr=  0.03,  y^»  0.6296,  and  y^=  0.1111. 

In  the  table  we  have  compared  plans  from  9  systems  defined  as  follows: 
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(1) .  Bayes.  Flans  minimizing  R(N,n,c). 

(2) .  IQL.  Min  R  for  P(pQ)  -  1/2,  pQ=  0.0250. 

(3) .  LTPD.  Min  R  for  P(p2>  -  0.10. 

(4) .  AQL.  Min  R  for  Q(px)  -  0.05. 

(5) .  Fixed  risk.  Q(p^)  *  0.05  and  P(p2)  *  0.10. 

(6) .  Percentage  inspection,  n  ■  0.057^+  0.1072»  0.04259. 

(7) .  Dodge.  The  AQL  system  with  5%  consumer's  risk  proposed  by  Dodge  in 

[6]  with  AQL  =  pj^. 

(8) .  Mil-Std.  Military-Standard  105D  with  AQL  =  p^. 

(9) .  Minimax  regret  plans  obtained  by  minimizing  max{K(N,n,c,p) -Nk  (p)’}. 

P 

For  each  of  7  lot  sizes  the  plans  and  the  corresponding  costs  have  been  found 
and  the  efficiency  has  been  computed. 

3.  Results  for  continuous  prior  distributions  of  process  average. 

If  the  process  average  varies  at  random  according  to  a  continuous  distribution 
with  density  w(p)  we  get  optimum  (Bayesian)  plana  with  properties  quite  different 
from  those  described  in  section  2.  The  asymptotic  results,  which  have  been  derived 
in  [4j,  are  the  following:  c  =  p^n  +  a  +  Ofa"1),  n  ■  +  X2+  o(l),  and  R  =  2n  +  o(l), 

where  a,X_,  and  \  depend  on  p  ,w(p  )  and  the  first  two  derivatives  of  w(p)  for 

p  *  pr* 

Furthermore,  the  asymptotic  efficiency  of  a  plan  (N,n^,Cj)  in  relation  to  the  optimum 
plan  (N,n0,cQ)  is 

•^.cp -/(£♦!*)  <i2> 

if  both  plans  use  the  right  relation  between  sample  size  and  acceptance  number, 
otherwise  the  efficiency  tends  to  zero. 

4,  The  minimax  regret  solution. 

If  the  prior  distribution  is  unknown  one  may  use  the  minimax  regret  method  to  derive 

a  sampling  plan,  see  [?] .  The  main  asymptotic  results  are  the  following; 

-1  2/3  1/3 

c  =  pfn  +  a  +  0(n  )  and  n  -  XjN  +  XjN  +  0(1),  where  a  depends  on  pf,  X^  and  X2 

on  p  and  the  cost  constants, 
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5.  A  more  general  model, 

i 

The  model  above  is  based  on  Che  assumptions  that  costs  are  linear  in  p,  and  that 
the  distribution  of  the  number  of  defectives  in  the  sample  (for  given  p)  is  binomial. 

In  the  present  section  we  shall  state  the  asymptotic  results  obtained  in  a  recent  paper 
[5j  concerning  a  much  more  general  model. 

Consider  a  problem  with  two  decisions  only,  acceptance  and  rejection  say,  where 

acceptance  is  preferred  if  the  unknown  parameter  0  is  less  than  or  equal  to  a 

specified  break-even  value  0 and  rejection  is  preferred  for  0  >  0 We  assume 

that  the  statistic  t  =  t(x  )  is  sufficient  for  0,  and  that  the  decision 

1  n  0 

rule  has  the  form:  "Accept  if  t  I  t  and  reject  otherwise."  The  problem  is  to  determine 
the  optimum  n  and  t^. 

We  shall  introduce  several  functions  of  0  such  as  w(0),  a(0),  etc.  and  to  shorten  the 
notation  we  write  w(0q)  =  w,  w ;(Cq)  “  w',  etc. 

The  model  is  based  on  the  following  assumptions: 

(1)  .  The  density  w(0)  of  the  prior  distribution  of  6  is  twice  differentiable  in 
an  open  interval  about  the  break-even  value  0^  and  w  >0. 

(2) .  The  loss  function  1.(0),  i.e.  the  loss  of  a  wrong  decision,  may  be  written 

as 

:'i1(8o-e)1(mu(a-0o)+  ...  )  toes  eQ 

1(0)  ..  ’  v  (13) 

I  I2(0-eo>  2(1+121(0-0O)+  ...  )  for  0>  e0 

where  arc  non-negative  and  (l^,^)  are  positive  constants. 

2 

(3) .  Conditional  on  0,  the  statistic  t  has  mean  0,  variance  a  (0)/n,  and 
coefficient  of  skewness  {3(0)/ i/n,  and  the  operating  characteristic  P(0)  of  the 
decision  rule,  i.e,  Pr(t  £  t®|0),  may  be  written  as  tax  Edgew*  rth  expansion 

P(0)  =  5>(x)-P(0)(x2-l)cp(x)/6/n  +  0(n_1),  x  -  (t°-0)Vn/a(0)  (14) 

where  cp  and  $  denote  the  density  and  the  cumulative  distribution  function  for  the 
standardized  normal  distribution.  The  functions  a(0)  and  {3(0)  are  assumed  to  be 
twice  differentiable. 

8  "1 

(4) .  The  (expected)  sampling  costs  are  r  (n)  *  kn  (l+0(n  )), 

s 

From  these  assumptions  we  find  the  average  regret  r(n,t®)  as 

r(n,t°)  =  rg(n)+r1(n,t°), 


(15) 


rv  m 
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where  the  average  decision  loss  is 


n  u 

r.(n,t°)  *  /  l(0)Q(0)w(0)d0  +  /  l(0)P(0)w(9)d0  , 


with  Q(0)  =  1-P(0). 

The  Bayesian  sampling  plan  (^q^q)  is  found  by  minimizing  r(n,t^).  Because  of  the 
continuity  of  1.(0),  P(0),  and  w(0)  in  the  neighbourhood  of  0^  and  the  decreasing 
standard  deviation  of  t,  the  asymptotic  solution  depends  only  on  the  properties  of 
the  functions  involved  at  0Q. 

It  is  obvious  that  we  might  have  chosen  one  of  the  costs  constants,  k  for  example, 
as  the  economic  unit,  i,e.  the  solution  will  depend  on  the  ratios  1^/k  and 
only.  In  deriving  the  asymptotic  results  we  assume  that  the  two  ratios  tend  to 
£  infinity,  and  that  l^/l^  is  constant.  We  shall  treat  n  and  t°  as  continuous  variables, 

From  (15)  it  follows  that  the  optimum  values  of  (n,t^)are  determined  from  the 

equations  dr./dt0  *  0  and  -dr./dn  -  dr  /dn.  For  any  n  the  optimum  value  of  t^,  i.e. 

—  0  —  9 

the  value  minimizing  r(n,t  ),  may  be  found  by  solving  the  first  of  these  equations. 
The  solution  will  be  denoted  tQ(n) .  Inserting  t°  ■  tQ(n)  int0  the  second  equation  we 

get  an  equation  in  n  alone  which  gives  r.  *  nQ  and  finally  tQ  “  to^no^  *  T*ie  func'cions 

ri(n,t0(n))  and  r(n,tQ(n))  will  be  denoted  r^n)  and  r(n),  respectively. 

Comparing  (10)  with  (15)  and  (16)  wc  find  that  (asymptotically)  the  model  simplifies 
to  the  one  considered  above,  if  the  cost  functions  are  chosen  as  r  (n)  ■  (p  -p  )n 

and  l{p)  “  N'|.P“Prl,  and  if  furthermore  a  (p)  *  pq  and  p(p)  «  (q-p)/i/pq,  since  (14) 

then  gives  the  Edgeworth  expansion  for  the  cumulative  binomial  distribution. 

It  is,  however,  clear  that  the  present  model  also  covers  cases  with  non-linear 
cost  functions,  and  operating  characteristics  different  from  the  binomial,  for 
example,  the  Poisson  and  the  normal. 


As  auxiliary  functions  we  introduce  the  incomplete  normal  moment  of  order  v 


m  (z)  *  /  tvcp(t+z)dt  -  /  (t-z)v  cp(t)dt, 

0  z 


v  §  0, 


Iv<z)  *  Ii  my(z)  +  i2  mv(‘z)- 

As  a  measure  of  efficiency  of  a  non-optimum  plan  (n, t^)  we  shall  use  the  ratio 


e(n,t°) 


r<V  t0)/r(n,t°). 
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This  measure  of  efficiency  lies  between  0  and  1,  it  is  invariant  to  the  choice  of 
origin  and  scale  of  costs  (or  utilities),  and  it  expresses  directly  what  costs  could 
be  reduced  to  (by  using  the  optimum  sampling  plan)  as  a  fraction  of  the  costs  for 
the  non-optimum  plan. 

If  t°  =. 0 j  +  o(l),  0^  4  ifc  is  easy  to  show  that  e  — >  0  as  n  — >  oo.  The  break¬ 
even  value  is  therefore  the  most  important  parameter. 


We  shall  in  particular  study  the  efficiency  first  under  the  assumption  that  tg(n) 
is  known,  so  that  the  efficiency  becomes  e(n)  =  r(ti0)/r(n),  and  next  under  the 
assumption  that  only  0Q  is  known,  so  that  t°  «  0Q  and  e(n,0Q)  ■  r(no,to)/r(n,0o)« 


In  stating  the  results  obtained  it  is  essential  to  distinguish  between  two  cases, 

v2  >  V1  say>  and  v2  “  vi  “  v* 


Results  for  v. 


v1  «  v. 


The  optimum  relationship  between  t  and  n  is  given  by 

t„  < 


0O  +  aa//n  4  bo/n  4*  0(n  ^^), 


where  a  and  b  are  determined  from  the  equations 


l.m  (a)  ®  l„m  (-a) 

-1  v  -2  v 


and 


b  =  -ocwVw  -  (v4l)a,4-(a  +v-l)p/6  -  pt  Uj jl  jin^  (a )+i2 il2mv+l  ^  “  a  ^ /-4+1  ^  ’ 

The  optimum  sample  size  is 

2/(2s4*v+l)  2s4*v4*l)/,„_ix  a  n/~"2/ (2s+v4-l), 


nQ-  7 


v4l 


14*287 


/  (2s4v4*l)  +  0(7 


•>) 


where  7  *  wot  ^^(a)/2sk  and 


6  »  a(ow'/w+(v4*l)p/6)4-  *2 


The  average  decision  loss  for  t  ®  tQ(n)  *s 


rx(n) 


;a 

v  4*  1  \'fnj 


v4l 


!  ,  +  (v+»g  +  0(I)  ] 
i1  +  (v+D.s  +  SV  /  ' 


and  the  minimum  regret  becomes 


,  v  ,  x  2s4*v+l 

r(n_)  ■  r  (n_)  — rr — 

'O'  sx  O'  v4*l 


1 

(*■ 


2s5 


\ 


(v4*2  )  (2s+v4*l )  vn0 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


■i 


m 

I 


For  the  efficiency  ve  get 


and 


1  „  v+1 

e(n)  2s+v+l 


+ 


2s  j  V| 

2s+v+l  \  n 


(v+l)/2 


1  v+1  ,n  .8 

e(n,0Q)'%/  2s+v+l  CnQ; 


+ 


28  /  !0  ; 

2s+v+l  n  ■ 


(v+1)/ 2 


W°> 
W>  ’ 


i.e,  the  second  member  of  (25)  has  been  changed  by  a  constant  factor. 


(25) 

(26) 


It  follows  that  tQ  — >  the  difference  tQ  -  0Q  being  of  order  1/ ,n  if  1^  =f  L, 

(non- symmetric  loss  function),  and  of  order  1/n  for  JL^  «  ^  because  then  a  =  0  from 

(20). 

If  ^  and  lg  are  proportional  to  lot  size  N  then  the  optimum  sample  size  n^  becomes 

proportional  to  jj2/(2s+v+1)  accor<j£ng  t0  (22).  By  varying  s  t  1  and  v  £  0  we  may 

1/2 

get  all  powers  of  N  less  than  or  equal  to  2/3.  For  v  ■  s  ■  1  we  get  N  . 

For  the  optimum  plan  the  ratio  of  the  average  decision  loss  to  the  sampling  costs 
tends  to  2s/ (v+1),  cf.  (2 :4). 

From  (25)  it  follows  that  the  efficiency  of  a  plan  based  on  the  right  relationship 
between  t®  and  n  will  be  high,  even  if  n  deviates  considerably  from  nQ.  For  v  »  s  *  1 
we  get:  (12). 

If  ©0  is  known  we  may  get  a  rather  high  efficiency  by  using  t°  ■  9^  even  if  n 
deviates  from  nQ,  see  (26). 

8  ? 

If  n  has  been  determined  from  a  wrongly  chosen  s  we  get  n  cc  V  8  4  °r  if  » 
discrete  prior  distribution  has  been  used  instead  of  the  correct  continuous  one 
we  get  n  a:  In  nQ,  and  in  both  cases  lim  e(n)  »  0. 

For  the  simple  symmetric  loss  function  .1(0)  ■  1^|0-—  6q|V  one  more  term  of  the  above 
expansions  may  be  found  in  [5|. 


Results  for  v2  >  v^. 

The  optimum  relationship  between  t°  and  n  is  given  by 

C0  "  +  ccm^'  “  oc((Vj+V2+l)ln  m  +  ln\  )/(m vTi)  +  o(l/m^  ft)  (27) 

where  m  ■  V(v2  -  v^ln  n  and  \  ■  (/2ff  a  1  1^ V(iir(vi  +  *))• 

The  optimum  sample  size  is 

2/  (2s+v  2+1 )  (v?+l )/  (2s+v?+l ) 

n0  -  7 2  ( (21n  72)/(2s+v2+1))  (1+o(1)),  (28) 

v2+1  (v2+1)/2 

where  y^  m  wot  (v2"V^)  l^2/2sk. 
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0 

The  average  decision  loss  for  t  =  tQ(n)  is 


v9+l  . 

r^(n)  ■  (wjL2/(^+l))(am/\/n)  L  (l+o(m’  )) 


and  Che  minimum  regret  becomes 


r(«o>  =  rg(n0)(2s  +  v2  +  l)/(v2  +  lKl+oCm^  )). 


For  the  efficiency  we  get 


.  v0*!-l  /  \S  „ 

1 _  ^  2  n  \  2s 

e(n)  *  28+v2+1  (n  Qj  +  2s+v2+l 


( noln  n  \ 
\  n  In  n„ / 


(v2+1)/2 


l  v2+1  / n_\s  2s  fn0ln  " 

!(n,eQ)  ~  28+Vj+I  \nQl  +  2a+v2+l  In  In  n( 


2s  !  n0ln  "  \(V2+1>/2  l(V2‘Vl>/2!  1lmv1+l(°)(V2+1) 


I  v9+l 


V  -V 

J  I2«  (Vj+1) 


so  that  lim  e(n,0  )  =  0  even  if  n  is  proportional  to  n 

u  u« 


For  v2  >  it  is  more  serious  to  reject  a  lot  which  should  have  been  accepted 
than  to  accept  a  lot  which  should  have  been  rejected.  (It  should  be  remembered  that 
for  large  n  it  is  only  losses  in  the  neighbourhood  of  0^  that  matter).  This  is  the 
reason  for  the  result  that  tQ  tends  to  0Q  from  above,  the  difference  tQ  -  0Q  being 
of  order  y  (In  n)/n,  so  that  the  rate  of  convergence  is  considerably  slower  than 
for  v2  -  v^. 

As  a  consequence  of  this  result  the  average  decision  loss  consists  essentially  of  the 

loss  from  wrongly  accepting  bad  lots.  For  the  optimum  plan  the  ratio  of  the  average 

decision  loss  to  the  sampling  costs  tends  to  2s/(v2+l),  see  (30). 
con- 

As  a  further  /sequence  it  will  be  seen  from  (32)  that  it  is  not  eatisfactory  to  put 
t^  ■  0q  because  such  a  plan  will  have  an  efficiency  tending  to  zero  for  n  — >  », 

The  results  above  depend  essentially  on  the  assumption  that  the  density  of  the 
prior  distribution  is  positive  and  differentiable  in  the  neighbourhood  of  0^,  and 
furthermore  on  the  assumption  that  1(0^)  ■  0  and  1(0)  >  0  for  0  4  0q. 

If  the  prior  distribution  is  discrete,  for  example  a  two-point  distribution  given 

1 

by  Pr(0®0. )  “  w. ,  i  =  1,2,  then  we  get  similar  results  as  in  section  2,  i.e. 
tQ  — >  0  where  0  depends  on  0^  and  ©2  only,  nQ  becomes  proportional  to  the 

logarithm  of  a  linear  combination  of  .1(0,)  and  1^09),  and  the  ratio  of  the  average 

~~  ^  ^  -1 
decision  loss  to  the  sampling  costs  will  be  of  order  n^  . 

If  1,(0)  «  0  for  0ji  0  ^  02,  i.e.  there  exists  an  indifference  zone  instead  of  a 
break-even  point,  then  we  get  similar  results  as  for  a  discrete  prior  distribution. 


Considering  restricted  Bayesian  plan3  it  is  clear  that  any  restriction  of  the 

form  P(t9.)  »  P,  say,  where  0.  4  and  P  is  a  given  number,  will  lead  to 
0  i  10 

t  — >  0^  and  therefore  to  an  efficiency  tending  to  zero.  Restrictions  should 
therefore  be  of  the  form  P(0Q)  =  P  or  PC^/QC^)  *  9>  where  p  is  a  constant  and 


®1  <  ®0  <  V 


6,  Discussion. 

Bayesian  sampling  plans  should  be  used  if  the  prior  distribution  is  known  and 
stable,  if  the  sampling  costs  and  the  decision  losses  are  Unown,  and  if  the 
purpose  is  to  minimize  the  sum  of  the  average  sampling  costs  and  decision  losses  for 
a  scries  of  lots. 

The  above  results  show  how  the  Bayesian  plans  depend  on  the  assumptions.  We  have 
found,  for  example,  how  the  optimum  sample  size  depends  on  the  main  parameters  in 
the  model.  Further  investigations  are  needed  to  evaluate  these  results.  It  is, 
however,  fortunate  that  the  formulas  for  the  efficiency  of  non-optimum  plans  show 
a  high  degree  of  insensitivity  to  deviations  from  the  optimum  sample  size 
if  only  the  break-even  point  is  known. 

In  practice  a  system  of  sampling  plans  is  often  required  to  serve  several  porposes 
-  besides  being  easy  to  administer.  In  particular  we  shall  here  mention  (a)  that 
the  system  should  protect  the  consumer  against  the  consequences  of  deterioration 
of  the  prior  distribution,  (b)  that  the  system  should  work  as  an  incentive  for 
the  producer  to  produce  better  quality  or  at  least  to  keep  to  the  quality  agreed 
upon,  and  (c)  that  average  costs  should  be  minimized.  So:  far,  however,  there  has  not 
been  developed  a  theory  taking  all  these  aspects  into  account. 

The  purposes  (a)  and  (b)  above  may  be  obtained  by  introducing  restrictions  on  the 
operating  characteristic,  and  by  alternating  between  normal  and  tightened  inspection. 
An  approximation  to  a  solution  may  therefore  be  obtained  by  using  a  restricted 
Bayesian  plan  with  decreasing  producer's  and  consumer's  risks  combined  with  a 
feed-back  mechanism  which  induces  shifts  between  normal  and  tightened  inspection 
according  to  changes  in  the  prior  distribution. 

It  would  be  very  useful  for  the  direction  of  further  research  in  this  area  if 
inspection  departments  would  publish  their  experiences  with  respect  to  costs  and 
prior  distributions.  Great  masses  of  data  must  exist  which  could  help  to  solve 
problems  regarding  the  form  and  stability  of  prior  distributions. 
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